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We consider the semiclassical ballistic cr-model as an effective theory describing the quantum me- 
chanics of classically chaotic systems. Specifically, we elaborate on close analogies to the recently 
developed semiclassical theory of quantum interference in chaotic systems and show how semiclassi- 
cal 'diagrams' involving near action degenerate sets of periodic orbits emerge in the field theoretical 
description. We further discuss how the universality phenomenon (i.e. the fact that individual 
chaotic systems behave according to the prescriptions of random matrix theory) can be understood 
from the perspective of the field theory. 
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I. INTRODUCTION 

In recent years, significant progress has been made in 
understanding the semiclassical basis of universality in 
quantum chaos. In a sequence of steps [TJ |5J |3J 2J [5J [5] , 
the theory of action correlations and periodic orbits has 
advanced to a stage, where the full information stored 
in the spectral correlation functions of random matrix 
theory (RMT) has become accessible by semiclassical 
methods. An important element in this development 
has been the establishment of a cross-link between semi- 
classics and the so-called zero dimensional sigma model 
(which in turn had long been known [7j to be in one-to- 
one correspondence to RMT): under conditions of global 
hyperbolicity and for low energies (energies lower than 
the inverse h/ts of the so-called Ehrenfest time, tg), in- 
dividual contributions to the semiclassical expansion of 
spectral correlation functions become fully universal in 
that they depend on combinatorics and topology of the 
underlying orbit correlations, but not on system specific 
details. Each such term has a corresponding contribu- 
tion to the loop expansion of the zero dimensional sigma 
model around one of its saddle points. From a certain 
perspective, the semiclassical analysis of spectral corre- 
lations has, thus, been tantamount to a reduction of the 
full microscopic information stored in the Green func- 
tions of an individual chaotic system down to the core 
data encapsulated in the zero dimensional sigma model. 

In this paper, we discuss an alternative reduction 
scheme, which is similar in spirit, but methodologically 
different. Our starting point will be the observation 
that the zero dimensional cr-model affords an alterna- 
tive interpretation, viz. as the globally uniform mean 
field limit of the ballistic cr-model. The ballistic cr-model 
in its variant considered here is an effective field theory 
of chaotic systems, obtained from a microscopic parent 
theory [5J EJ EH] by projection onto the sector of fields 
fluctuating on length scales > H. We will present vari- 
ous arguments to the effect that this projection captures 
the essential physics of the problem. The resulting effec- 
tive field theory is defined on shells of conserved energy 
in classical phase space. Its dynamics is controlled by 
the - equally classical - Liouville operator. Quantum 



mechanics enters the theory through a feature known 
as 'non-commutativity'. In practice, non-commutativity 
implies that points in phase space can be defined only up 
to a maximal resolution set by the Planck cell. (Notice 
that the typical linear extensions of Planck cells are of 
0{h 1 / 2 ), much bigger than the cutoff length of the the- 
ory.) Below we will explore the conspiracy of classical 
hyperbolicity and Planck cell resolution in the long time 
dynamics of chaotic systems. 

Specifically, we will show that a perturbative expan- 
sion of the model leads to structures that are remark- 
ably similar, if not equivalent, to those arising in peri- 
odic orbit theory: taking the first quantum correction 
to the spectral correlation function (the 'Sieber-Richter 
term' [2]) as an example, we will recover information on 
the phase space available to individual families of pe- 
riodic orbits, and the corresponding action correlations. 
The actual integrals describing these corrections turn out 
to be identical to those of periodic orbit theory. This find- 
ing suggests a quantitative identification of the Feynman 
diagrams of the theory with the action-correlated orbit 
pairs of semiclassical analysis (However, an extension of 
the calculation to higher orders in perturbation theory, 
necessary to substantiate this claim, is beyond the scope 
of the present paper.) 

Going beyond perturbation theory, we will demon- 
strate that for correlation energies smaller than the in- 
verse of the Ehrenfest time non-uniform fluctuations in 
phase space are effectively suppressed. The resulting the- 
ory of uniform fluctuations (viz. the zero dimensional 
cr-model) predicts universal behavior, in agreement with 
the predictions of RMT. Field theory thus provides a 
comparatively simple way of understanding the basis of 
the Bohigas-Gianonni-Schmid universality hypoth- 
esis, alternative to the explicit summation over infinite 
hierarchies of periodic orbits [B] . The computational ef- 
ficiency of the 'mean field plus fluctuations' approach to 
universality arguably represents the most important ad- 
vantage of the present formalism. 

The rest of the paper is organized as follows: in section 
|n]we review the semiclassical ballistic cr-model. In sec- 
tions |lll] and |IV] repsectively, we explore the field theory 
approach to universality, and the connections to periodic 
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orbit theory. A number of issues relating to the deriva- 
tion of the theory will be addressed in appendix [A] 

II. BALLISTIC (t-MODEL 

We wish to explore correlations in the level density, 
p(E), of globally hyperbolic (chaotic) quantum systems, 
as characterized by the two point correlation function 
R 2 (s) = A*((p(E+^) - (p(E+^)))(p(E) - 0,(25)))). 
Here, (...) denotes averaging over E over an interval 
[E - E av /2,E Q + E av /2], where E av < Eq. The goal 
is to show that for energies uj — sA/n < Ht^ the func- 
tion i?2(s) approaches the spectral correlation functions 
of RMT. Here, t E = A -1 ln(c 2 /^) is the Ehrenfest time, 



Eq. (J2j) has the status of an effective field theory in classi- 
cal phase space, obtained from a microscopic parent the- 
ory - the energy averaged field integral representations 
of the microscopic Green functions [U [10] - by elimi- 
nation of field configurations fluctuating on short scales 
< h. In Appendix |A| we will argue that Eq. ^ relates 
to the microscopic formalism in much the same way the 
semiclassical Gutzwiller sum relates to the microscopic 
Feynman path integral. 

The integration variables in T(x) = {T QQ '(x)} 
are matrix valued fields defined on shells T = {x|_ff(x) = 
Eq} of constant energy in classical phase space. Here, 
x = (q, p) where q and p are coordinates and momenta, 
respectively, H (x) is the Hamiltonian function of the sys- 
tem, and the integral over the energy-shell is normalized 
to unity, J r (dx) — 1. For time reversal and spin rotation 
invariant systems (orthogonal symmetry class, (5 = 1), 
the 'internal' structure of the matrices T aa is described 
by a composite index a = (a,r,t), where a = +/— dis- 
criminates between the advanced and retarded sector of 
the theory, r = 1, . . . , R is a replica index, and t = 1, 2 
accounts for the operation of time reversal. Time rever- 
sal symmetry reflects in the relation a^ 1 ' T T = T -1 , 
where er^ are Pauli matrices and the superscript 'tr' in- 
dicates action in time reversal space. For time rever- 
sal non-invariant systems (unitary symmetry, /3 — 2) no 
time-reversal structure exists and a = (a,r). In either 
case, the matrices T carry a coset space structure in 
the sense that configurations T and TK are identified 
if [K, A] = 0, where A — erf 1 ' and 'ar' stands for action in 
advanced/retarded space. 

The fluctuation behavior of the fields T in pi is gov- 
erned by the classical Liouville operator {H, } (where 
{ , } is the Poisson bracket.) Quantum mechanics enters 



A the dominant Lyapunov exponent of the system and c 2 
some classical reference scale of dimension 'action' whose 
specific value will not be of much relevance. [T2] 

We represent the spectral correlation function in terms 
of a replica generating functional (Choosing the replica 
variant of the theory is motivated by its high suitabil- 
ity to perturbative calculations [3]; it is a matter of a 
straightforward redefinition of the field target space to 
upgrade the formalism to a supersymmetric field theory), 

R 2 (s)=~Be Km-^d 2 s Z(s), (1) 

where R is the number of replicas and the generating 
function is defined as 



(2) 
I 

the problem through the back-door, viz. by the presence 
of Moyal products in (|2j) . The Moyal product between 
two phase space functions A and B is defined as 

P) 

where / = ( °j jj ^ is the symplectic unit matrix. (For all 

practical purposes, the definition ^ will be more conve- 
nient than the standard representation [TS] , (A * B) (x) = 
exp(f <9J,/<9 X ) A(x')B(x)\ x=x ,.) The presence of the 
Moyal product implies that (a) all quantities appear- 
ing in the theory get effectively averaged (smoothened) 
over Planck cells. Rclatedly, (b) the generator of classi- 
cal time evolution {H, } acts on distributions smooth in 
phase space on scales ~ W , rather than on mathematical 
points. In the following sections we will explore how the 
interplay of hyperbolic dynamics {{H , }) and Planck cell 
smearing (*) determines the output of the theory. 



III. UNIVERSAL LIMIT 

We wish to explore the behavior of the theory ([2| for 
correlation energies u ~ 2-Kht^ 1 fit^ 1 of the order 
of the inverse Heisenberg time tn = 2irhA~ 1 and much 
smaller than the inverse of the Ehrenfest time. For sim- 
plicity, we will consider systems with broken time reversal 
invariance throughout this section. 

For energies u 3> A s 3> 1, the partition func- 
tion Z(s) may be evaluated by perturbative methods. 
To prepare the perturbative expansion of the action, we 



Z(s) = J DT exp(-5[T]), S[T] — ^(dx)tr^T -1 * A{H. T} + AT -1 * AT 
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introduce the rational parameterization 



T = t + W, W 



-fit 



B 



(4) 



where the block structure is in advanced/retarded space 
and the generators B are Rx R complex matrices. Sub- 
stitution of this representation then leads to a series 
S[B, Pt] = V~ =1 S^ 2 ") [B, St] t where is of nth order 
in B and P^. Specifically, 



S<®[B,B*] =t H J{dx)tr{B^L u B), 



where 



-iu/h~{H, } 



(5) 



(6) 



and the Moyal product stars have been omitted for muta- 
tional simplicity. Higher order terms in the action con- 
tain traces over (Moyal) products of matrices BB' . . .BK 
Due to the Planck-cell 'averaging' inherent to the product 
([3]) , the Wick contraction of individual matrix elements 
of B and B^ will generate expressions 



/>oo 

P,(x,x') = / ^e iw * fi " I P t (x ) x / ) 
Jo 



where 



P., 



Pi 



J{H, } 



(7) 



(8) 



are the (retarded) Liouville propagators in the energy and 
time representation, respectively, x and x' are arbitrary 
points in phase space and 



1 

w 



^/(x + y) 



is symbolic notation for coordinate averaging over a 
Planck cell. We interpret P t (-,x') as the dynamical 
evolution of a smooth phase space distribution of exten- 
sion fv and centered around x'. It has been rigorously 
shown by methods of symbolic dynamics |14j that for 
times t < £e> these distributions are centered around the 
classical trajectory through x', i.e. the dynamics is ap- 
proximately described by the Liouville evolution of x'. 
However, beyond t ~ t-g,, the distribution rapidly (over 
time scales comparable to t mlx ) crosses over to a uniform 
distribution in phase space. These structures are sum- 
marized in the ansatz 



P t (x,x')^{^ (x 



where = is the volume of the energy-shell, and the 
normalization is such that for our unit-normalized phase 
space integral, f r {dx) Pt(x, x') = 1. 

For energy scales oj < ht^} the dominant contribu- 
tion to the time integral in d7| comes from large times 



*'(*)), 



t<t E , 
t > t E , 



(9) 



t E < t < hui^ 1 . One expects that in this regime phase 
space fluctuations will have effectively relaxed to spatial 
uniformity. To describe the damping of inhomogcneous 
modes in more explicit terms we follow a prescription for- 
mulated by Kravtsov and Mirlin [15 , then in the context 
of diffusive systems: employing the ansatz 



T(x) =T (l + W P (x)) 



(10) 



we isolate the inhomogeneous contents of the fields T. 
Here, To = const, describes the zero mode sector and 

Hp(x) = W(x) - Wq(k), Wq(x) = J (dx) W(x), 

is a projection onto purely fluctuating field configura- 
tions: J(dx)Wp(n) — 0. Substitution of Eq. (10 1 into 
the action and expansion to second order in Wp — 

^ B ^ generates the decomposition 

S[T , W P ] = 5 [Qo] + S&[B p , B P ] + 5 C [T , B P , Bp], 

(11) 

where 



2A 



tr(Q A), Q = r AP- 1 (12) 



is the zero mode action, and 



SW[B p ,B P ] =t H J(dx)ix(B\,CoB P ), 



(13) 



the quadratic action of the fluctuating fields governed by 
the generator £ — £w=o = S — {H, }, S \ of the time 
integrated dynamics. Finally, 



S c [Tq, Bp, Bp] 



2h 



(dx) x 



(14) 



x tr ((T - 1 AT ) 11 P P P p - (T - 1 AP ) 22 P P P P ) (x), 

where the superscripts refer to advanced-retarded space. 
Following Ref. 15, we expand in S Cl retaining only con- 
tributions of minimal order in to: 



+{S c [T ,Bp,Bp)) - ^((S c [T Q ,B P ,Bp}) 2 ), 

where (...) = / D(B p , Bp) exp(-M 2 ) [B p , B P ]) (. . . ). 
While the contraction of the P's in the first contribu- 
tion to the second line vanishes in the replica limit, the 
second term generates the effective action 



S eS [Q ] = S a [Q ] - ^MQqA)) 2 x 

x |(dx)(dx / )(Pp)o(x,x / )(Pp)o(x' ! x), (15) 
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where, again, the subscript P stands for projection onto 
the fluctuating sector: 

(F P )o(x,x') = (l - J (dx)) (l-J (dx'))p (x,x') 

di (<5(x-x'(i)) -fT 1 ) , (16) 

and in the last line we have switched to a time represen- 



tation. Notice that the Liouville propagator in Eq. (15) 
is evaluated on single phase space points, x and x', rather 
than on Planck-cells x and x'. This is because the phase 
space integral of the Moyal-product of two operators col- 
lapses to the ordinary product [16], i.e. 

J (dx) tr ((T^ATo) 11 ^ * B p ) (x) 

(dx)tr ((T - 1 AT ) 11 4(x)B p (x) 



Substituting the result Eq. (16) into Eq. (151, we obtain 



S e s[Qo] — Sq[Qq] 



16H 2 



(wtr(QoA)) 2 x 



f dt /(dx) (S(x-x(t)) -ft' 1 ) . (17) 
Jo Jr 



According to Eq. (17 1, the correction term is given by the 



integrated weight of periodic orbits x — ► x of duration t, 
minus the total phase volume. Using that [T7] 

lim [ dt [ (dx) S(x - x(t)) = T, 

one concludes that this term cancels, i.e. to lowest or- 
der in the expansion, fluctuations of the inhomogcneous 
modes do not change the universal zero-mode action. At 
higher orders in the expansion, we are met with expres- 
sions (cf. Eq. d7l|) 



(PpUx,x') = 
= £ dte^ (l-J (dx)) (l-J (dx'))P t (x,x') 



tE 



die' 



T(W( x -x'(i))-l), 



where Eq. (|9| was used. The important point here is 
the truncation of the time integral at tE- This means 
that at nth order in the expansion in fluctuation prop- 
agators, Pp, corrections ~ (orfe/ft)™ will be generated. 
Our analysis above shows that terms with n < 3 cancel, 
i.e. corrections to the zero mode action can arise only at 

o^tEh- 1 ) ~ o(t E /t H ) 3 , 



the leading correction to the universal result scales as 

(ut E /hf. 

Summarizing, our analysis shows that at correlation 
energies uj < H/t-E inhomogeneous fluctuations get effec- 
tively damped out and the field theory collapses to an in- 
tegral over the zero mode. This means that RMT results 
(plus weak corrections in the parameter (wtE / 'h) n - 3 ) will 
be obtained for spectral correlation functions, and other 
observables probing the long time behavior of the system. 
Our 'mode damping' approach to universality is comple- 
mentary to semiclassical analysis 6 , where the spectral 
correlation functions are constructed explicitly, by sum- 
mation over infinite hierarchies of periodic orbits. The 
connection between these approaches will be explored in 
the next section. 



IV. PERTURB ATIVE EQUIVALENCE TO 
SEMICLASSICS 

Having shown that the ballistic er-model crosses over to 
the universal zero dimensional model at low energies, we 
now approach the problem from a different perspective. 
We will perform a straightforward perturbative expan- 
sion around the high energy saddle point A to elaborate 
on parallels between the field theory ^ and the peri- 
odic orbit approach to spectral correlations. Focusing 
on the lowest order quantum corrections to the spectral 
form factor (the so-called Sieber-Richter term [2]), we will 
show that the two formalism are structurally similar, to 
an extent that the correlations of individual periodic or- 
bits can be reproduced from the field theory formalism. 

Throughout this section, we will consider time reversal 
invariant systems, i.e. (3=1. 



A. Semiclassics 

For the benefit of non-expert readers, we begin with a 
brief review of recent semiclassical results. Consider the 
Fourier transform of the spectral correlation function, 



KM = - / dsR 2 (s)e 

7T 



-2isr 



(18) 



The RMT result for the (3 — 1 form factor reads as 
[ 2t-t1h(1 + 2t), t < 1 



T > 1. 



(19) 



Specifically, the short time expansion of the spectral form 
factor in t < 1 (corresponding to the expansion of the 
spectral correlation function in ui S> A) starts as 



S cS [Qo] = S [Qo] + 0(wi B /fi)"- 3 



K{t) = 2r - 2r 2 



(20) 



This observation is consistent with the periodic orbit 
analysis of Ref. [TB] where it has been shown that 



In semiclassics one aims to reconstruct that expansion 
from the Gutzwiller double sum over periodic orbit pairs 
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(a, a'), 



K{t) 



eUs«-s a >) d l T 



2t 



H 



(21) 

where (...) is an average over orbit energies (and a small 
interval of orbit times.) and A a , S a , and t a are the sta- 
bility amplitude, the action, and the time of traversal of 
orbit a, respectively. 



The first term in (20 1 obtains from Berry's diagonal 
approximation , 

K^\t) = 2r, 

i.e. an approximation that retains only identical or- 
bits a = a' and mutually time reversed a — a' orbits, 
and uses the Hannay-Ozorio de Almeida []15] sum rule 
( E Q l^a| 2 ^( r — ta/tn)) — t to determine the weight of 
the remaining orbit sum. 





t E 



FIG. 2: Phase space representation of the encounter region. 
Discussion, see text. 



Imagine a Poincare section cutting through the en- 
counter region where the participating orbit segments 
have their (avoided) self intersection (cf. Fig. [2]) Choos- 
ing one of the piercing points, x, as the origin of a lo- 
cal coordinate system (cf. inset of Fig. [2|, we intro- 
duce 2/ — 1 = 3 energy shell coordinates (t, u, s), where 
t is a time-like coordinate running along the trajectory 
through x, and u/s are coordinates along the locally un- 
stable/stable directions (the existence of the latter guar- 
anteed by the assumption of global hyperbolicity.) One 
may then show [51 |5U] that the contribution of all orbit 
pairs containing a single encounter region is given by 




FIG. 1: Cartoon of an orbit pair contributing to the first 
quantum correction to the spectral form factor. 

The structure of the leading order (in r) quantum cor- 
rections to the spectral form factor was first identified 
by Aleiner and Larkin pQ: two orbits differing in the 
presence or absence of a small angle self intersection in 
configuration space interfere to provide a stable contribu- 
tion to the double sum (cf. Fig. [I]) For the sake of later 
comparison to the tr-model formalism, we briefly review 
the quantitative computation of this contribution [2] in 
the invariant language of hyperbolic phase space coordi- 
nates [20 for a system with / = 2 degrees of freedom. 



duds t 



2t enc (u, s) 



(22) 

Here, tenets) = A 1 ln(c 2 /ws) is the time of traver- 
sal through the encounter region, A the (self-averaging) 
Lyapunov exponent of the system, and c 2 some classical 
action scale. We note that (i) the dominant contribu- 
tions to the integral are of order us ~ h implying that 
t cnc (u, s) ~ A^ 1 ln(c 2 /h) is of the order of the Ehrcnfest 
time, (ii) the product us is an invariant of the motion 
(Liouville's theorem), which means that the choice of the 
Poincare section entering the construction is arbitrary, 
an d (ii i) the only surviving contribution to the integrand 
in (22) are of order 0(t enc (u, s)°). (Terms containing 
t cnc , either from the numerator, or the expansion of the 
denominator can be shown to effectively oscillate to zero 
in the semiclassical limit [3j|4j[6].) 

In later work, the perturbative analysis of action cor- 
relations has been extended first to cubic and then to 
infinite order [31 U |S] in the r-expansion. Employing an 
unconventional representation of the two-point correla- 
tion function [6] it has also been possible to reconstruct 
the bottom part of Eq. (19 1, which is not accessible in 



terms of a straightforward series expansion of the con- 
ventional two-Green function representation. 



6 



B. Field Theory 



Eq. (22) encapsulates the phase volume available to 



individual pairs of interfering trajectories (pre-exponcnt) 
and the corresponding action correlation (exponent). Ex- 
tending the analysis of Ref. 16, we will show below that 
the same information is stored in the effective action ^ . 
Although our discussion will be restricted to first order 
in perturbation theory in r (the extension of the anal- 
ysis to higher orders in perturbation theory is left for 
future work), it suggests a general equivalence between 
periodic orbit theory and the perturbative expansion of 
the cr-model. 

In the time reversal invariant case considered here, the 



blocks B in Q are (2R x 2P)-matrices subject to the 
time reversal symmetry condition B* = — <7 2 r B T a 2 ■ 

Substitution of this representation into the action ^ 
generates the first order correction to the spectral corre- 
lation function 



{d*)(tY{B^BB^CuB)) 



where we have omitted the Moyal product stars for the 
sake of notational simplicity, and averaging is over the 
quadratic action, (...) = J V{B\ B) e -s <2) [s,s t ](. .. ). 
Doing the Gaussian integrals over matrix elements of B, 
we obtain (cf. Ref. [TH] for technical details), 



n 2 

Re — d 2 s 



dyidy 2 

(27rfi) 2 (/- 1 ) 



:* y = 'yip w ( x + yi/2,x - yi /2)£^ x+y2/2 P w (x + y 2 /2, x - y 2 /2). 



Here, x = (q, — p) T is the time reversed of the phase space point x = (q, p) T , the integrals over yi and y 2 represent 
the generalization of the Moyal product ^ to the trace of a product of four operators and the subscript x in £^ iX 
indicates on which coordinate of the propagator P w (x,x') the derivatives of the Liouville operator act. To make 
further progress, we Fourier transform the expression above whereupon it transmutes to the first quantum correction 
to the spectral form factor, 



-2r 2 D 2 / (dx) 



dyidy 2 e^ Iyi 
(2nh) 2 U-i) 



dtfPt-t>(* + yiAx - yi/2)A,, x+y2/2 P t ,(x + y 2 /2,x - y 2 /2),(23) 



r 



Here, C t = d t — {H, } and P t is the propagator in time 
representation (cf. Eq. pj.) We proceed by introduc- 
ing a coordinate system that has x as its origin and 
x+x' <-* (r, u, s), where r is the coordinate of x' along the 
classical trajectory running through x and u and s pa- 
rameterize the components of x' in the locally stable and 
unstable direction around x (see Fig. [2]) We assume [2T] 
that the function P t (x + y/2,x — y/2) exhibits the fol- 
lowing characteristics: (i) Pt(x + y/2,x — y/2) = Pt(u) 
depends only on the unstable component of y. Indeed, 
this component controls the rate at which the trajectory 
starting at x — y/2 deviates from that through x. This 
deviating component (rather than the approaching com- 
ponent, s) determines the spatial and temporal structure 
of the Lyapunov region. Conversely, P t (x+y/2,x — y/2) 
depends only on s. (ii) In order for Pj(x + y/2, x — y/2) 
to become non-vanishing, the two trajectories through 
x — y/2 and x + y/2 must have left the Lyapunov region 
around x (cf. the 'large scale picture' Fig. [TJ. This takes 
a time of order 2 x t(u), where 

t{u) = A _1 ln(c/u) 



accounts for 'half of the encounter time. Additional to 
this time, some time of classical duration, roughly of or- 
der i m i x , is required to 'tie' the outgoing and incoming 
trajectory segments to a closed link. These assumptions 
are summarized in the ansatz 



P t (x + y/2,x-y/2) =n~ 1 Q(t -2t(u)), (24) 



where is a smeared step function interpolating between 
zero and unity over a time interval of order t m - lx . Nei- 
ther the detailed structure of O, nor the exact value of 
t(u) are of further relevance. We note, however, that 
the postulated independence of P t of the 'longitudinal' 
coordinate, r, implies stationarity of the long time prob- 
ability distribution Pt>t(u) under the Hamiltonian flow, 
{H,P t>t{u) } = 0. 



Using that C t Pt 
becomes (/ = 2) 



dtPt, the first quantum correction 
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Ki2) ( T ) = -J^§j2 J dt ' J du l ds 1 du 2 ds 2 e^^- s ^P t . tl {s l )d t) P tl {u 2 ) 
02 L dt ' J dudseiUSPt - t '^ d ^) Pt '^ 



T 2 n 

2ttH 

2jrh, 
2r 2 f 

— / dudse* us d tcnc{u , s) 

^ [dudse^- 2 *™^. 
2nh J 2t enc (u,s) 



r 



dudse* us d t ( u )(t - 2t{u) - 2t(s)) 
t - 2t cnc (u, s) 



(25) 



in agreement with the semiclassically derived expression 
(22 1. In the first line we used that, under the condi- 



tions stated above, the long time probability distribu- 
tion is invariant under the action of the Liouville op- 
erator. In the fourth line, we introduced the full en- 
counter time, t cnc = t(u) + t(s), and the last equality is 
based on the above mentioned vanishing of any power 
(t enc (u, s)) n upon integration against the oscillatory ker- 
nel ~ exp(iH~ us). 




FIG. 3: Assignment of propagator coordinates (lines termi- 
nating in thick points) to Feynman amplitudes (lines contin- 
ued through the Poincare section). 

The derivation above demonstrates that field theory 
and semiclassics, resp., allocate the same phase space 
volume to single encounter processes. We finally show 



an interpretation in terms of individual periodic orbits. 
To this end, we consider a Poincare section through the 
encounter, as shown schematically in Fig. [3] The four 
lines terminating in the two pairs of connected points in 
the phase space plane represent segments of classical tra- 
jectories beginning at x + y 12 /2 and x — y lj2 /2, resp., 
i.e. idealized classical trajectories corresponding to the 
arguments of the field theory propagators in ( 23 1 . Each 



that the field theoretical expression ( 25 1 indeed affords 



of these points has a stable (s) and an unstable (it) co- 
ordinate. The fact that the propagators P are retarded 
implies that we can distinguish between 'incoming' and 
'outgoing' terminal points. Now, consider the behavior 
of a trajectory carrying the stable coordinate of an in- 
coming terminal propagator point, and the unstable co- 
ordinate of a terminal point of the other propagator pair 
(cf. inset of Fig. [3]). The trajectory running through 
this reference point will interpolate between the propa- 
gator stretches involved, and is naturally interpreted as 
part of a closed loop. The analogous definition of three 
more points (see inset of Fig.[3| plus associated trajectory 
stretches leads to the identification of a periodic orbit 
and its topologically distinct partner orbit. The action 
difference between these two orbits is but the product 
us. This construction implies an interpretation of the 
field theory expression ( 25 1 purely in terms of periodic 
orbits. The correspondence between the two formalisms 
is established before integration over phase space coor- 
dinates (u,s,t): apparently, the action ^ encapsulates 
detailed information about the hyperbolic dynamics and 
action correlations of individual Feynman amplitudes. 



V. SUMMARY AND DISCUSSION 

In this paper we have considered the semiclassical bal- 
listic <7-model as an effective theory of quantum chaos. 
Defined in classical phase space, the semiclassical a- 
model is a field theory whose dynamics is driven by the 
classical Liouville operator. Quantum mechanics enters 
through a structure known as non-commutativity. In 
practice, non-commutativity (i) limits the maximal res- 
olution of the theory to structures of the order of the 
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Planck cell, and (ii) leads to the appearance of charac- 
teristic phases, which play a role analogous to the action 
correlations of periodic orbit theory. 

We have shown in perturbation theory that the a- 
model describes spectral correlations in far reaching anal- 
ogy to semiclassics. The fact that semiclassics and field 
theory attribute the same weight to individual orbit cor- 
relations in phase space (cf. Eqs. (22 1 and (25 1) makes 



one suspect that the a- model fully encapsulates the infor- 
mation carried by the Gutzwiller double sum. However, 
this expectation has not been proven beyond first order 
in perturbation theory. 

Perhaps most importantly, we have shown that the 
semiclassical cr-model provides for an efficient description 
of the crossover to the universal regime of RMT correla- 
tions at energies below the inverse of the Ehrenfest time. 
Adapting a technique originally developed by Kravtsov 
and Mirlin to describe spectral correlations in disordered 
systems, we saw that in chaotic systems — in contrast 
to disordered systems — quadratic inhomogeneous fluctu- 
ations do not give relevant corrections to the universal 
action. Extending this approach to higher orders in the 
inhomogeneous modes, a systematic quantitative esti- 
mate for the corrections to RMT spectral statistics in 
individual chaotic systems can, in principle, be obtained. 
Methodologically, this way of approaching the universal- 
ity phenomenon is complementary to (and arguably more 
economical than) the infinite order summations over cor- 
related orbit pairs of Refs. |3j |U 

We have enjoyed fruitful discussions with P. Brouwer 
and F. Haake and thank S. Heusler for critical 
reading of the manuscript. Work supported by the 
Sonderforschungsbereich SFB/TR12 of the Deutsche 
Forschungsgemeinschaft . 



APPENDIX A: REGULARIZATION 

The cr-model discussed in the foregoing sections has the 
status of an effective theory, obtained from an underlying 



'bare' theory by elimination of rapid fluctuations. In this 
section we discuss the status of that reduction. Com- 
paring to semiclassics, we will argue that the projection 
onto an effective field theory resembles the passage from 
the Feynman path integral to the semiclassical Gutzwiller 
sum. 

The action of the bare ballistic cr-model [5] |5J HO] is 
given by 



S q u[T] 



2£ , tw 



tr ( T~ 1 A[H, T] 



As 
2tt 



at -i at 



(Al) 

where, H is the Hamiltonian operator of the system and 
tr is a trace over Hilbert space, projected to an energy 
strip of with T av centered around Eq- The integration 
variables T = {T"3 } are operators in a product Hilbert 
space spanned by real space coordinates, q, and 'internal' 
coordinates, a. (The internal structure of the matrices 
T qq ' has been discussed in section [Ii] above.) 

The action Sq U [T] is obtained from the rigorous func- 
tional integral representation of the two level correla- 
tion function after (a) a saddle point approximation 
(which generates the nonlinear constraint Q 2 = 1, where 
Q = TAT -1 ) and (b) first order expansion incommu- 



A2 



tators [H,T 1 ]. As we will argue in section 
these two approximations arc largely immaterial. 



below, 



Semiclassical Model 



with the (semi) classical con- 
made by switching to a 
Upon Wigner transforma- 



ble 



Preliminary contact 
tents of the theory is 
Wigner representation, 
tion of Hilbert space operators A c 

J dAqe lH lAq ' p Aj +Aq / 2 , q _Aq/2, the operators -> 
T(x) transmute to the fields in classical phase space dis- 
cussed above. |H] Further, tr(. . . ) — > J r g?^x (...) be- 
comes an integral over a phase space energy shell ^E av , 
centered around Eq and of thickness T av . We are thus 
led to the phase space representation [16] 



SAT] = 



i/3ir 
2T a v 



tr IT" 1 *A[H * T] - - AT" 



AT 



(A2) 



r 



where the asterisks stand for Moyal products, as usual. 



Importantly, the field theory (A2) does not contain 



mechanisms inhibiting the buildup of rapid fluctuations. 
The action S[T) of fields fluctuating on quantum scales 
~ h a , a > 1 in classical phase space is qualitatively dif- 
ferent from the semiclassical action considered above: for 
fields fluctuating transverse to the energy layer of thick- 
ness -E av ~ fr/to ~ h reduction of the phase space integral 



to an integral over a single energy 'shell' is no longer pos- 
sible. Further, the series expansion of the Moyal product 

[H*T] = ih{H, T} + h 3 0(d 3 Hd 3 T). 

shows that for fields fluctuating on scales ~ h a , a > 1, 
the approximation of the quantum commutator by the 
Poisson bracket is no longer permissible. 
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A save way to eliminate those rapid field fluctuations is 
by adding a weak random potential to the Hamiltonian. 
An ensemble average over that randomness will generate 
a second order differential operator which may be tailored 
so as to effectively remove fast field fluctuations. As dis- 
cussed in Ref. [23] a 'quantum random potential' V q ~ h, 
vanishing in the classical limit (and parametrically in h 
weaker weaker than the 'regulator' suggested by Alcincr 
and Larkin [TJ suffices to remove field fluctuations 
on phase space scales ~ h. 

Less rigorously, one may argue that the action S of 
rapidly fluctuating fields will lead to highly oscillatory 
exponents exp(iS) which likely tend to average out: a 
glance at Eq. ^ shows that the largest contributions 
to the action of the problem, set by the largest value 
of the correlation energy w ~ h/t m i x , scale as A" 1 x h. 
However, fields fluctuating on scales ~ h a , a > 1 will lead 
to contributions of A -1 x h n ^ a \ where the odd index 
n designates the order of the Moyal expansion. One may 
argue that in the semiclassical limit, these terms generate 
strong phase cancellations which render the contribution 
of strongly fluctuating fields effectively meaningless. (For 
a caveat in this argumentation, see section A 2 below.) 



We note that the above phase cancellation argument 
resembles the logics inherent to the stationary phase 
derivation of the Gutzwiller trace formula from the Feyn- 
man path integral; there, too, avoidance of rapid fluc- 
tuations is the dominant principle. More specifically, 
Gutzwiller's trace formula is based on a stationary phase 
approximation to the Feynman path integral in S/h 
(where S stands for the typical value of a classical ac- 
tion.) This approximation effectively averages over fine 
structures on scales h a , a > 1. Relatedly, the semiclas- 
sical analysis of spectral correlations will be oblivious 
to the averaging of the system Hamiltonian over 'quan- 
tum' perturbations V q ~ h. (Not affecting the classical 
dynamics, such perturbations merely change the phases 
weighing individual periodic orbits. In the evaluation of 
the Gutzwiller double sum of spectral correlations, these 
phases cancel out. [25] ) 

Summarizing, the semiclassical cr-model (|2]) considered 
in the main text obtains as a projection of the 'quan- 
tum' cr-model (Af I onto the sector of fields fluctuating 
on scales > H. This projection may be effected either by 
averaging the system over a quantum random potential 
of strength ~ h, or by alluding to the prospected irrele- 
vance of strongly fluctuating actions in the semiclassical 
limit (i.e. ad hoc restriction of the functional integral to 
fields fluctuating on scales > H.) 

One thus reduces the microscopic theory (equivalent to 
the full Feynman path integral representation of Green 
function correlators) to an effective theory (likely equiv- 
alent to the Gutzwiller approximation of the Feynman 
path integral.) 

2. Zero modes 

While the above phase cancellation arguments apply 
to 'generic' field configurations, there is one family of 



rapidly fluctuating configurations that deserves separate 
consideration: the quantum action (Al) is nullified by 
a large number of zero modes |26] fluctuating at scales 
of the order of the Fermi wave length. Choosing a rep- 
resentation in terms of cigenfunctions, H\ip a ) — e a \tp a }, 
wherein T = {T$} and [H,T] = (e Q - e' a ){T^}, we 

conclude that there exist N ~ E av /A zero modes 
whose action is controlled only by the energy contri- 
bution, S qn [T aa ] = ^trCAT-^AT^). Upon Wigner 
transformation, the zero mode operators turn into N zero 



mode functions of the action (|A2j , rapidly fluctuating in 
classical phase space. These modes are certainly not 'un- 
physical'. In the case of time reversal non-invariant sys- 
tems, one may indeed verify |27j that the integration over 
T ao leads to the formally exact eigenvalue decomposition 



Ra(u) 



J2 6{u-{e a -e b )). (A3) 



|£o-e a ,6l<-E av /2 



What makes the rigorous derivation of Eq. ( A3 1 possible 



is a mathematical principle known as 'semiclassical ex- 
actness'. [25] (Incidentally, we note that the possibility 
to obtain an exact representation of the spectral correla- 
tion function from the quantum ballistic cr-model shows 
that the approximations on which the derivation of the 
latter is based - saddle point approximation and first 
order expansion in the quantum commutator - become 
exact in the limit -E av /A — > oo.) 



In spite of its formal correctness, Eq. ( A3 1 is, of course 



useless in practice (much like the evaluation of the Feyn- 
man path integral in an exact basis of eigenfunctions, for- 
mally an exact alternative to the semiclassical stationary 
phase approximation, would be pointless.) 

The conservative way to remove the zero modes is, 
again, by averaging over weak randomness. Indeed, semi- 
classical exactness represents a very delicate structure; 
even minisculc changes in the action will spoil the exact 
cancellation of all non-Gaussian fluctuations on which 
Eq. (A3 1 is based. (Although we are lacking a rigorous 



justification, we believe that averaging over a random po- 
tential whose inverse scattering time, or level broadening, 
is as weak as fir^ 1 ~ A will suffice to effectively remove 
the zero modes.) ft is, then, favorable to switch to a 
Wigner phase space representation, and classify fluctua- 
tions along the lines of the semiclassical analysis of sec- 
tion |A l| Alternatively one may avoid averaging over ran- 
domness and accept the presence of zero modes - after all 
the integration over these modes produces a meaningful, 
if useless result. Mapping the integral onto a phase space 
representation restricted to slowly fluctuating modes one 
may deliberately sacrifice the exact information stored in 
the zero modes in return for a semiclassically meaningful 
approximation scheme. 
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